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Birth, Death and Flight: A Theory of Malthusian Flocks 
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I study "Malthusian Flocks": moving aggregates of self-propelled entities (e.g., organisms, cy- 
toskeletal actin, microtubules in mitotic spindles) that reproduce and die. Long-ranged order (i.e., 
the existence of a non-zero average velocity {v(f,t)) ^ 0) is possible in these systems, even in spa- 
tial dimension d — 2. Their spatiotemporal scaling structure can be determined exactly in d = 2; 
furthermore, they lack both the longitudinal sound waves and the giant number fluctuations found 
in immortal flocks. Number fluctuations are very persistent, and propagate along the direction of 
flock motion, but at a different speed. 

PACS numbers: 05.65.+b, 64.70.qj, 87.18.Gh 



Flocking [l[ - the coherent motion of large numbers of 
organisms - spans a wide range of length scales: from 
kilometers (herds of wildebeest) to microns (microorgan- 
isms @, 0] ^mobile macromolecules in living cells (3, HI)- 
It is also [6[ a dynamical version of ferromagnetic or- 
dering. A"hydrodynamic" theory of flocking [7| shows 
that, unlike equilibrium ferromagnets Q, flocks can spon- 
taneously break a continuous symmetry (rotation in- 
variance) by developing long-ranged order, (i.e., a non- 
zero average velocity (v(f,t)) ^ 0) in spatial dimensions 
d = 2, even with only short ranged interactions. 

Many quantitative predictions of the hydrodynamic 
theory , including the stability of long-ranged order in 
d = 2, the existence of propagating, dispersion-less sound 
modes with non-trivial direction dependence of their 
speeds, and the presence of anomalously large number 
fluctuations, agree with numerical simulations]?! H| , and 
experiments on self-propelled molecules (To|. 

However, a recent re-analysis ll| of this hydrodynamic 
theory has cast doubt on the claim that exact scaling ex- 
ponents could be determined for flocks in d = 2. This is 
due to the erroneous neglect inQ of non-linearities aris- 
ing from the local number density dependences of various 
phenomcnological parameters; these non-lincaritics could 
change the scaling exponents from those claimed in Q ■ 

In this paper, I show that these difficulties can be 
avoidedin flocks without number conservation. A num- 
ber of real systems [l^, [l^ lack number conservation, in- 
cluding growing bacteria colonics [13(, and "trcadmilling" 
molecular motor propelled biological macromolecules in 
a variety of intra-cellular structures, including the cy- 
toskleton, and mitotic spindles [1, in which molecules 
are being created and destroyed as they move. Hence, he 
study of such systems is not only convenient, but experi- 
mentally relevant. The most obvious example of flocking 
- namely, actual birds - is clearly not a good example 
of a Malthusian flock. I will henceforth use the term 
"boid" [l[ to refer to the moving, self-propelled entities 
that make up the flock. 

I will treat systems with the same symmetries as were 
considered in earlier work on immortal flocksQ: orien- 
tationally ordered, translationally disordered phases (i.e., 
phases with (v(r,t)) ^ that are uniform in space and 



time) in systems with short-ranged, rotation invariant 
interactions, moving through (or on, in d = 2) a fixed 
background medium that breaks Galilean invariance. 

Since my treatment is hydrodynamic, it only describes 
large systems at long length and time scales. However, 
it becomes asymptotically exact in that limit. 

The removal of number conservation leads to profound 
changes. The sound modes of number conserving (here- 
after "immortal") flocks disappear, and are replaced by 
longitudinal velocity fluctuations which drift in the di- 
rection of flock motion with a speed 7 7^ vq, where vo is 
the mean speed of the flock. While drifting, these modes 
also spread diffusively along the direction of flock motion, 
and hyperdiffusively perpendicular to that direction. 

In both Malthusian and immortal flocks, anomalous 
hydrodynamics stabilizes long-ranged orientational order 
(i.e., (v(r, t)) ^ 0) in spatial dimension d = 2. In Malthu- 
sian flocks, the order outlives the boids: the persistence 
time diverges as the number of boids N — s- 00, while the 
lifetime of the boids remains finite in this limit. 

The scaling exponents of the hydrodynamics can be de- 
termined exactly in spatial dimension d — 2 for Malthu- 
sian flocks, while, as discussed above, those in immortal 
flocks cannot. These exponents are: the dynamical ex- 
ponent z for the scaling of timescales t{L±) oc L\ with 
length scale L± perpendicular to the direction of flock 
motion; an anisotropy exponent £ for the scaling of dis- 
tances L\\(L±) oc parallel to the direction of motion 
with L± ; and a "roughness" exponent \ relating the 
scale of velocity fluctuations to L± via Sv oc . I find 
that, for Malthusian flocks in spatial dimension d = 2, 
these exponents arc: 
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(1) 



The velocity field can have long-ranged order; 
((v(f,t)) ^ 0) in d = 2, because the roughness exponent 
X(d = 2) <0. 

Number fluctuations in Malthusian flocks exhibit 
anomalous persistence: the experimentally observable 
density-density correlation function 



C p (f, t) = (Sp(f, t')Sp(r' + r, f + t)> 



(2) 
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where Sp(r, t) = p(r, t) — po is the departure of the local 
number density of boids p{r,t) from its mean value po, 
decays algebraically with time at a fixed point in space: 



C p {f= 0,i) oc |i| 



(3) 



in spatial dimensions d = 2, while for a point translat- 
ing along the direction xii of flock motion at the "drift" 
speed 7 the decay is even slower: in both 2 and 3 spatial 
dimensions, I find 



C p (r = ftx\\ , t) oc \t\~ 



(4) 



These should be contrasted with the exponential decay 
with time of density fluctuations that occurs in a disor- 
dered Malthusian flock (i.e., one in which < v >= 0). 

Since C p and C v can be constructed from any set of 
high-resolution images of a moving flock, as has been 
done for flocks of starlings in reference fl4ll and for a re- 
lated correlation function of bacteria in[13||. it should not 
be difficult to test these predictions. 

I will now outline the derivation of these results. 

My starting equation of motion for the velocity is ex- 
actly that of an "immortal" flock[l5[: 



d t v + X 1 (v-V)v + A 2 (V • v)v + A 3 V(|??| 2 ) = 

av- P\v\ 2 v- VPi -v(v- VP 2 (p, \v\)j 

+D° B V(W ■ v) + D T V 2 v + D 2 (v-V) 2 v + f (5) 

where all of the parameters \{i = 1 — > 3), a, /3,D B , Dt2 
and the "pressures" Pi, 2 (/9, \v\) are, in general, functions 
of the boid number density p and the magnitude \v\ of 
the local velocity. I will expand Pi 2 (p, \v\) about po- 
P(p) = PP + J2Zi cn,n(\v\)Sp n , where i = 1, 2. 

In ©J (3, D B , D 2 and Dt are all positive, while a < 
in the disordered phase and a > in the ordered state. 

The a and (3 terms give v a nonzero magnitude vq = 

in the ordered phase. The diffusion constants Db,t,2 
reflect the tendency of "boids" to follow their neighbors. 
The / term is a random Gaussian white noise, mimicking 
errors made by the boids, with correlations: 

< Mr, t)fj(P, t) >= ASijS^r- ?)8{t - t') (6) 

where A = constant, and i,j label vector components. 
The "anisotropic pressure" P 2 {p, \v\) in §5§ is only al- 
lowed due to the non-equilibrium nature of the flock; in 



an equilibrium fluid such a term is forbidden by Pascal's 
Law. In earlier work [?J this term was ignored. 

Note that ([5]) is not Galilean invariant; it holds only in 
the frame of the fixed medium through or on which the 
creatures move. 

I now need an equation of motion for p. In immor- 
tal flocks, this is just the usual continuity equation of 
compressible fluid dynamics. For Malthusian flocks, it 
must also include the effects of birth and death. I will 
assume that the death rate goes up as the density goes 
up (the "Malthusian" assumption [16j ) . so that the dif- 
ference between the two - that is, the net, local growth 
rate of number density in the absence of motion, which 
I'll call g(p) - vanishes at some fixed point density po, 
with larger densities decreasing (i.e., g{p > po) < 0), and 
smaller densities increasing (i.e., g(p < po) > (J). 

The equation of motion for the density is now simply: 



dtp + V • (vp) = g{p) 



(7) 



Note that in the absence of birth and death, g(p) = 0, and 
equation 10 reduces to the usual continuity equation, as 
it should, since "boid number" is then conserved. 

Since birth and death quickly restore the fixed point 
density po, I will write p(r, t) = po + 6p(r, t) and expand 
both sides of equation Q to leading order in 6p. This 
gives po\/ -v = g'(po)Sp, where I've dropped the d t p term 
relative to the g'(po)5p term since I'm interested in the 
hydrodynamic limit, in which the fields evolve extremely 
slowly. This equation can be readily solved to give 



AD) 



g'(po) 



01,1 



-(V-z?) 



(8) 



where AD^- 1 is a positive constant, and eri.i is the first 
expansion coefficient for Pi . I can now insert this solution 
(0) for Sp in terms of v into the isotropic pressure Pi; the 
resulting equation of motion for v is: 

d t v + X 1 (v-V)v + A 2 (V • v)v + A 3 V(|t7| 2 ) = 
av- [3\v\ 2 v- v Off- VP 2 (p, \v\)j 

+P4 1) V(V • v) + D T V 2 v + D 2 (v-V) 2 v + f, (9) 

where I've defined = D° B + AD$ . Taking the dot 
product of both sides of ([9]) with v itself, and defining 
U(\v\) =a{\v\,p)- (3(\v\,p)\v\ 2 , I obtain: 



- (a t |«| 2 + (Ai + 2\ 3 ){v- V)|it| 2 ) + A 2 (V • v)\v\ 2 = U(\v\)\v\ 2 - \v\ 2 v ■ VP 2 + D^v- V(V • v) 

+ D T v-V 2 v + D 2 v-({v-V) 2 vj+v-f, (10) 
I 

In the ordered state (i.e., in which (v(r,t)) = vqX\\), and I can expand the v equation of motion for small 
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departures 5v(f, t) of v(f, t) from uniform motion with 
speed v : 

v(f,t) = (v + Sv ll )x l] +v x (r,t) , (11) 

where, henceforth || and _!_ denote components along and 
perpendicular to the mean velocity, respectively. 

In this hydrodynamic approach, I'm interested only in 
fluctuations Sv(r,t) = 6v\\x\\ + v±(f,t) and 8p(f,t) that 
vary slowly in space and time. Hence, terms involving 
spatiotcmporal derivatives of 8v(r,t) and Sp(r,t) are al- 
ways negligible, in the hydrodynamic limit, compared to 
terms involving the same number of powers of fields with 
fewer spatiotemporal derivatives. Furthermore, the fluc- 
tuations 8v(r, t) and Sp(f, t) can themselves be shown to 
be small in the long-wavelength limit. Hence, we need 
only keep terms in |10| up to linear order in 5v(f, t) and 
5p(f, t). The v ■ f term can likewise be dropped. 

These observations can be used to eliminate many 
terms in equation [10] . and solve for the quantity U = 
(a(p, \v\) - P{p, \v\)\v\ 2 )- I obtain: U = A 2 V • v + v ■ VP 2 . 
Inserting this expression for U back into equation [5], I 
find that P 2 and A 2 cancel out of the v equation of mo- 
tion, leaving, ignoring irrelevant terms: 

d t v + Xi(v- V)t7+ A 3 V(|t/| 2 ) = D T V 2 v 
+D { g ] V(V -v) + D 2 {v-W) 2 v + f, (12) 

This can be made into an equation of motion for v_i_ in- 
volving only v± (r , t) itself by projecting perpendicular 
to the direction of mean flock motion x\\, and eliminat- 
ing 5v\\ using U = A 2 V ■ v + v ■ VP 2 and the expansion 

U « -r 1( 5v|| - r 2 <5p, where I've defined T x = - {J^j 

and r 2 = — (j§f^ , with super- or sub-scripts de- 
noting functions of p and \v\ evaluated at p — po and 
\v\ = Vq. Doing this, and using ((5J) for p, I obtain: 

d t v± + jd\\v± + Xi(v± ■ V±)v± + A 3 V_l (|w_i| 2 ) = 
D T \7 2 ± v ± + D B V ± (V ± ■ v ± ) + D l{ dfv ± + f u (13) 

where I've defined 7 = \\v$, Db = D B + 2uoA3(A 2 — 

r 2 AD^Vo-i)/ri and D\\ = D% + D%v 2 . 

Changing co-ordinates to a new Galilean frame f ' mov- 
ing with respect to our original frame in the direction of 
mean flock motion at speed 7 - i.e., r' = r — ^tx\\ - gives 

d t v± + \i(v ± ■ V ± )v ± = D T V 2 ± v ± + D s V_l(V± • v±) 

+Pi,a^7 ± + /i . (14) 

Ignoring the non-linear term Ai in this equation of mo- 
tion gives a noisy, anisotropic, vectorial diffusion equa- 
tion. This can be readily solved for the mode structure 
and fluctuations by spatiotcmporal Fourier transforma- 
tion, and has d — 1 diffusing modes in spatial dimen- 
sion d. These separate into d — 2 "transverse" modes 



(i.e., modes with v± perpendicular to q±, all with the 
same imaginary eigenfrequency: — icu T = Dx\q ± \ +D\\q 2 . 
The remaining diffusive mode (the only mode in d = 2) 
is "longitudinal" (i.e., has v ± along q ± ), with frequency 
-iu L = DjifJ 2 + D\\q 2 , where D ± = D B + D T . 

Because the dynamics described above is in the 
Galilcanly boosted frame, the dynamics in the original 
reference frame r will have a steady drift at velocity 7 
superposed on the diffusive motion described above; that 
is, both eigenfrequencies get jqu added to them. 

I can also calculate the real-space velocity fluctuations 
{\v±(r,t)\ 2 ) in this linearized approximation; I find that, 
in this approximation, this diverges in all d < 2. This 
is analogous to the Mermin- Wagner theorem Q in equi- 
librium magnets. However, as in immortal flocksQ, this 
"Mermin- Wagner" result, and all the linearized scaling 
laws, arc invalidated for d < 4 by the Ai term in (TT4")) . 

To show this here, I'll analyze equati on (1141) using the 
dynamical Renormalization Group(RG)|17|. 

The dynamical RG starts by averaging the equations 
of motion over the short- wavelength fluctuations: i.e., 
those with support in the "shell" of Fourier space 6 _1 A < 
\q\ < A, where A is an "ultra-violet cutoff', and b is an 
arbitrary rescaling factor. Then, one rescales lengths, 
time, and v ± in equation (|14p according to v ± = b x v' L , 
r± = br[_, rj| = 6^(r'|)', and t = b z t' to restore the ultra- 
violet cutoff to A. This leads to a new equation of motion 
of the same form as (|14[) , but with "renormalizcd" values 
(denoted by primes below) of the parameters given by: 

D' BtT = b z - 2 (D B . T + graphs) , (15) 
D\ =6 Z - 2 <0D|| + graphs) , (16) 
A' = b z -i- 2 x +1 - d (A + graphs) , (17) 

A'i, 3 = & z+x - 1 (Ai,3+ graphs) , (18) 

where "graphs" denotes contributions from integrating 
out the short wavelength degrees of freedom. If we ig- 
nore these graphical corrections (valid for Ai^ small), and 
choose z, Cj and x to keep the linear parameters D By T,\\ 
and A fixed, equation (|18[) implies that an initially small 
Ai will grow for all d < 4, meaning the linearized theory 
It is possible to get exact exponents in d = 2. This 
is because the nonlinearitics - the Ai and A3 terms - in 
(fnj add up to a total derivative in d = 2 (specifically, 
+ A3) d±v\), since the _L subspace is one dimensional 
in d = 2. In contrast, in immortal flocks, v-p non- 
linearities arising from the p-dependence of Ai cannot be 
written as total derivatives, making it impossible to ob- 
tain exact exponents, a fact missed by Q. Here, because 
the Ai term is a total _L-derivative, it can only graph- 
ically renormalize terms involving _L-derivatives them- 
selves. Hence, the graphical corrections to D\\ and A in 
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equations (|16j) and (jTTJ) vanish. Hence, at a fixed point, 
in d = 2, 

z - 2C = 0, z - C - 2 X + 1 - d = z - C - 2 X - 1 = 0. 

(19) 

There are no graphical corrections Ai either, because 
the equation of motion (|14p remains unchanged by the 
transformation: f± — > f j_ — Xiilit , vj_ — > v± + v\ 
for arbitrary constant vector V\ _L xm. This exact sym- 
metry must continue to hold upon renormalization, with 
the same value of Ai. Hence, Ai cannot be graphically 
renormalized. Requiring that X[ = Ai in (j 18[) . and set- 
ting graphs = 0, implies \ = 1 — z in all d < 4. This 
and (|19|> forms three independent equations for the three 
unknowns x, z, and £, whose solution in d = 2 is ([1]). 

The scaling exponents z, £, and x determine the scaling 
form of the velocity-velocity autocorrelation function in 
arbitrary dimension d through the scaling relation^: 



C v (r , t) = (v ± (0, 0) • v ± (r , i)) = |^ ( ^!L_ , * 



\r , | 2x G 



r-jj - 7* t 



r, \* 



(20) 



where the second equality follows from scaling arguments 
applied to the boosted equation of motion (|14[) . and the 
third arises from undoing the boost. Here G(u, w) is a 



-7< 



and 



scaling function, with scaling arguments it = ,:, 
w = jp^ji- The asymptotic limits of G(u, u>) and C v {r, t) 
can be obtained by the following arguments. 

When ru — -ft — > and f — 5- 0, C v (f,t) must clearly 
depend only on rj_, and should not vanish. Hence G(u <C 
<C 1) — > constant 7^ 0. This in turn implies that 
C v (f,t) oc I 1 2x for |rj c > |r|| — 7*] , t*. Similarly, 
if r — >• and t — > 0, then C v (r,t) should depend only 
on |?*|i — jt\. This implies G(u,w) <x u < for u ^> w, 1, 



in order to cancel off the |^| 2x prcfactor in Equation 

(|2^1) . This in turn implies that C v (f,t) oc |rii — 7t|T~ for 

|rii — jt\ 3> , Similar reasoning implies that 

C v (f,t) oc |t|^ for f 3> Ir^l 2 , |rii — jt\i . Hence, using 
the exact exponents ([!]) in d = 2, 



C„ (f, t) oc 




|rj.|s > |r|| -7t| , 7:2 
|t-|| - 7t| > , ^ (21) 
1*1 > |rj.|* , hi -7*| 2 ■ 



This correlation function can be measured directly in 
both simulations 0, E|, and experiments [l4^ . 

The relation © between density and velocity implies 
that density correlations should obey the same sort of 
scaling law, but with an additional power of for 
every power of Sp; hence, in d = 2: 



C p (r,t) 



'II _ 7* 



t 




rj_p |rj_| 
\r±\i > |r|| - 7*| , 
hi - 7*1 > |n_ I* , ti 



(22) 



The last line holds in d — 3 as well, because x = 1 ~ z j 
does. The last two lines of (|2"2l directly imply equations 
(|3]) and Q . It can also be shown [ll| that at equal times 
C p (f,t = 0)decays sufficiently rapidly that there are no 
giant number fluctuations in Malthusian flocks. 
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